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Although gauge-boson propagators in asymptotically-free gauge theories satisfy a dispersion re- 
lation, they do not satisfy the Kallen-Lehmann (K-L) representation because the spectral function 
changes sign. We argue that this is a simple consequence of asymptotic freedom. On the basis of the 
QED-like Ward identities of the pinch technique (PT) we claim that the product of the coupling 
and the scalar part d{q^) of the PT propagator, which is both gauge-invariant and renormalization- 
group invariant, can be factored into the product of the running charge <f'{(f) and a term H{(f) both 
of which satisfy the K-L representation although their product does not. We show that this behavior 
is consistent with some simple analytic models that mimic the gauge-invariant PT Schwinger-Dyson 
equations (SDE), provided that the dynamic gauge-boson mass is sufficiently large. The PT SDEs 
do not depend directly on the PT propagator through d but only through H. 

PACS numbers: 11.15.Tk, 11.15.Kc UCLA/09/TEP/42 

I. INTRODUCTION 
A. Positivity and gluon mass 

A long-standing problem of non-Abelian gauge theories (NAGTs) is the lack of positivity of the imaginary part 
of the gauge-boson propagator, violating the Kallen-Lehmann (K-L) representation. This was first pointed out in 
an early paper on the gauge- invariant pinch-technique (PT) propagator Later, many authors found the same 
behavior in lattice simulations of the gauge-dependent and unphysical propagator of the Landau gauge 0, 0] (and 
references therein). Although it is questionable to assign a physical meaning to this lack of positivity in a gauge- 
dependent quantity such as the Landau-gauge propagator, many authors see it as a sign of confinement, since the 
propagator of an unconfined field presumably has a normal K-L representation. Aubin and Ogilvie Q] trace it to 
technical deficiencies in lattice gauge-fixing procedures. 

In this paper we argue that it is plausible (but unproven) that this lack of positivity is an elementary consequence 
of asymptotic freedom, and is simply resolved by a factorization of the PT propagator into two terms each of which 
is both gauge- and renormalization group (RG)-invariant, and each satisfies the K-L representation. We construct 
some simple analytic models of the PT propagator and vertex that illustrate the necessary positivity and absence 
of unphysical singularities, provided that there is a sufficiently-large dynamical gluon mass. We also discuss models 
of the Analytic Perturbation Theory (APT) type, which can satisfy positivity with zero gluon mass but still have 
unphysical behavior. 

Section fll Al briefly covers notation as well as some background on the PT. In Section fll Bl we argue that the lack 
of positivity is a simple consequence of asymptotic freedom and the fact Q that the product of the coupling 
and the (scalar part of the) gauge-invariant PT propagator d{q^) is not only gauge-invariant but also RG-invariant, 
independent of the choice of a renormalization point. This makes this product a truly physical quantity. The same 
is true for the photon propagator in QED, as has been known for decades, and for the same reason: The Ward 
identities of QED or of the PT require that the gluon vertex function renormalization constant and a wave-function 
renormalization constant be the same. 

Not every non-perturbative approximation for the propagator can be expected to satisfy these positivity constraints, 
in an asymptotically-free gauge theory. In the rest of the paper we construct non-perturbative models that do satisfy 
them, provided that there is a large enough dynamical gluon mass. 

Section III CI is an illustration, within the context of an analytically-soluble model similar to an earlier [l| one- 
dressed-loop Schwinger-Dyson equation (SDE) for the PT propagator, of how the positivity argument above can only 
be realized with a sufficiently large gluon mass. In the model, provided that the gluon mass m is large enough (on 
the QCD scale A) each factor in the product (f'H behaves precisely as would be expected, with no bizarre behavior 
coming from non-positivity. But if m/A is less than a critical value mc/A spurious singularities arise, such as ghost 
or spacelike poles in the propagator. We estimate mc/A ~ 1.2 in our study of the one-dressed- loop PT propagator; 
given the approximations made there, we believe the range should be from 1 to 1.5 or so. This effectively provides. 
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as we will see, an upper limit to the running charge at zero momentum: OsfO) = ^^(0)/47r < 0.5 — 0.7. This is fairly 
consistent with other determinations from phenomenology 3, studies [J, [a, Sill of the PT SDE, and a study of the 
functional Schrodinger equation Q. 

Section IIIII briefly reviews the evolution of APT from an originally massless form [10| with correct positivity 
properties, yet showing unphysical behavior, to a massive form rather similar to that of [l| and the present paper. In the 
original APT positivity was satisfied, although the PT upper limit on as{0) is exceeded, with as{0) ~ l/(47r6) « 1.1. 
This in itself is not necessarily serious, but what is serious is that the APT running charge, even though finite at 
zero momentum, has infinite slope. A later work 11] corrects this deficiency by invoking an ad hoc gluon mass, in 
somewhat the same spirit as the PT gluon mass, with results quite similar to our first model. Other authors [12] have 
also invoked masses as cutoffs for APT. 

Section IIVI remarks on the important fact that the vertex SDE can be reformulated entirely in terms of the 
propagator factor H with a positive spectral function and a special half-proper vertex that is both gauge- and RG- 
invariant; the original PT propagator d, with its positivity violations, never appears. This reformulation avoids possible 
violations of positivity that could allow unphysical vertex behavior. We illustrate with an analytic approximation 
inspired by a one-dressed-loop toy model [l^ of the Schwinger-Dyson equation for the three-gluon vertex and show 
that it has a spurious spacelike singularity if the gluon mass is too small. 

Section HV CI is a discussion of certain typical all-order extensions and resummations of the massless toy model that 
still lead to unphysical singularities; these can only be resolved with a dynamical gluon mass. 



B. The positivity problem and asymptotic freedom 

Write the obvious factorization of the product g^d: 

gH[q^)^f{q^)H[q^) (1) 

where g{q^) is the gauge-, scheme-, and renormalization-point-independent running charge of the PT. The other factor 
H{q^) has the same properties, since the product does. We argue that both factors obey a standard K-L representation 
with a positive imaginary part (our metric is such that q^ > for timelike vectors): 

fie) = r da^^ (2) 



Jim^ (y -q^ -le 

and p, pH are positive. The lower limit involves the dynamical gluon mass m, which we discuss later. 

We plausibly know the behavior at infinite momentum of both factors in the product. Asymptotic freedom tells us 
that 

^('^'^ 61n(4vA2)- 

where b is the lowest-order coefficient in the beta-function and A is the QCD scale. As for H{q^) there is no reason 
from perturbation theory or non-perturbative PT constructions to believe that it departs from the simple free-field 
behavior for large momentum. For example, it is well-known that the PT is equivalent order by order to the 
background field Feynman gauge 14 1, and old perturbative calculations in this gauge through two loops [l^ show 



that all large-momentum logarithms are accounted for in the running charge. The result is that at large momentum 
the PT propagator vanishes according to d{q^) ~ l/(g^lng^). But a propagator vanishing more rapidly than 1/q^ 
implies that the spectral function in the would-be K-L representation is necessarily negative somewhere. There is, 
we claim, only an indirect connection — at least for the gauge-invariant PT propagator — between non-positivity and 
confinement (a connection only to the extent that asymptotic freedom implies confinement). In fact, the gluon is not 
confined, but screened, in the usual sense that the string in the adjoint-representation Wilson loop always breaks at 
sufficiently large distance even with no adjoint matter fields. We will not attempt any analysis of non-positivity in 
the Landau-gauge, but it is likely that even if some non-positivity comes from the Aubin-Ogilvie 3] effect, there will 
still be some residual non-positivity coming from asymptotic freedom. 

What happens at infrared momenta? It has long been argued [ll,[i^ that the QCD gluon should pick up a dynamical 
mass that completely preserves local gauge symmetry. This is consistent with phenomenology (for example, (sl. fl^). 
and a number of studies of the PT Schwinger-Dyson equations [Ts^ have found P, Hll] a PT pole mass m of order 
0.6 GeV. These studies also indicate that the mass runs with momentum and should be denoted m{q^), consistent 
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with the operator-product expansion result that m?{q^) vanishes (modulo logarithms) like (G uuG^^) I at large 
momentum [l9l . There is a large body of lattice-simulation evidence in the Landau gauge [20n2ll. [23 . [23l . [2^ . [25l . 
m, [13, H [SliO, [3l|, [si, [H [M, [Si, [si, [53, [si for a gluon mass of several hundred MeV. (The pole mass of the 
gluon propagator in any gauge is gauge-invariant and physical, although because it is a timelike pole it is not easy to 
determine from lattices simulations.) We think it plausible, then, that H{q^) has some such form as: 

H{q^) = 2 (4) 

m^(q ) — q-' - le 

where m{q^) is the running mass. For simplicity and brevity, we use in this paper a fixed gluon mass; running does 
not interfere with the main positivity arguments. Then _ff is a simple free massive propagator. 

II. A MODEL OF THE PT INVERSE PROPAGATOR 
A. A few words on the pinch technique 

Begin with some notation. The pinch-technique propagator has the form: 

A„0(g)=P„^((7)rf(g)+e^; PM=-9c.P + ^- (5) 
The corresponding inverse pinch-technique propagator is: 

Kli^) = PMil^ +"^{1)] + -^q^qp. (6) 

The scalar function d is completely independent of the gauge chosen. 

The PT is a systematic way of extracting gauge-invariant proper self-energies, vertices, and the like from gauge- 
invariant quantities such as the S-matrix. The PT propagator, for example, is not constructed just from the usual 
Feynman graphs; it also receives contributions from other graphs through a so-called "pinch" , in which longitudinal 
momenta in numerators, coming from vertices and propagators, trigger Ward identities leading to the replacement of 
certain propagators by unity. This changes the topology of the graphs where this occurs, and some of these pinch 
contributions are readily recognized as contributions of propagator type. Although not recognized at the beginning, 
it is now known [ij] that the PT is the same graphical expansion as that of the background field-Feynman gauge. 
Because the whole point of the PT is to maintain gauge invariance, it is essential that, when this graphical expansion 
is resummed to a dressed-loop or skeleton expansion, all Green's functions appearing in the skeleton expansion obey 
the correct Ward identities. In the case of the PT these are the naive Ward identities of QED, with no ghosts. 

It would appear that progress can only be made by solving all possible SDEs at once, since that will guarantee 
satisfaction of the Ward identities. However, it is possible to find approximate three- and higher-point proper vertices 
that satisfy the Ward identities exactly and that are expressed solely in terms of the PT propagator itself. This 
approximation, known as the gauge technique, is valid for infrared-dominated phenomena. Although we will not give 
any details here, we have in mind the gauge technique of Ref. which gives the following expression (group indices 
suppressed) for the gauge technique proper vertex: 

where the first term on the right is the bare vertex F*^. Here H^j, = P^^^H is the PT proper self-energy introduced 
above. This vertex satisfies: 

fcr{f „/3^(fcl, ^2, /eg) - f°^^(fci, ^2, fca)} = tli3^{k2) " 6/3^(^3) (8) 

no matter what the choice of H is. In consequence, an approximate but fully gauge-invariant one-dressed-loop SDE 
for the PT propagator can be written solely in terms of that propagator. It is in this sense that we speak of studying 
the PT propagator on its own terms, without further specification of the vertex beyond that of the gauge technique. 
In this paper we use the self-contained PT propagator equation of Ref. [i|. Later we will study a toy model of a 
"half-proper" three-gluon PT vertex SDE in which some factors coming from the propagators are incorporated in the 
definition of the vertex, and we arrive at essentially self-contained vertex SDEs. 
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B. Sign problems 

Define a K-L function as a real-analytic function with at most one pole and a cut along the real positive axis, 
satisfying an unsubtracted dispersion relation with a positive spectral function. 

We conjecture that the PT propagator is the product of two K-L functions. The product of two K-L functions may 
or may not be a K-L function, but in the present case we know it cannot be, because the product vanishes faster than 
1 1 near infinity. The product of two K-L functions having this property therefore requires a spectral function that 
changes sign. Multiply together two K-L functions, call them G\ and G^- For each function we have: 

G.= / do^^ (* = 1,2) (9) 

with Pi nowhere negative. The product G1G2 obeys the dispersion relation: 

G,G, = / daP-l^ (10) 
J q^ -a 

with 

PiM = v [ + (11) 

J a - a' 

which may be negative in places. 

We note parenthetically that it is possible, in certain field theories involving scalar particles, for the propagator 
(taken to be K-L) to have a single zero between the particle mass and the lowest two-particle threshhold and for 
the proper vertex to have a pole at the same place. In our case this would correspond to a zero in H{q^) for 
< q^ < Am?. But for propagator models what we use as a criterion for a critical mass, based on asymptotic 
freedom, is not to exclude a zero of the propagator but to exclude an unwanted pole. It is true that the vertex models 
we study define a "critical" mass by excluding vertex singularities, but these are not related in any obvious way to 
zeroes of H . In any case we explicitly exclude the possibility of a zero in H by assumption, and our techniques show 
no signs of such a zero developing. See [s^ for further details and references. 



C. Mass parametrization 

The conventional approach to the propagator, whether from the PT or elsewhere, is to calculate in some approx- 
imation the proper self-energy, that is, the inverse propagator. The technique of going from the propagator to the 
inverse propagator (or vice versa) reminds us of analytic perturbation theory (APT), which we discuss in the follow- 
ing section. APT was used in the NAGT context to render 01^(0) finite, even in the zero-mass limit. However, an 
unphysical singularity remains, and in fact mass-improved APT also has a critical mass that is rather close to the 
values we give in this section. 

We assume that there is at most one pole in the PT propagator, at q^ = m?, representing dynamical gluon mass 
formation, and no zeroes, and we replace the factor H of Eq. (g]) by the simple massive propagator (g^ — -t- ie)~^, 
where the mass does not run. Our pretensions to accuracy in the infrared do not justify saving the running, so m? 
can be thought of as either the running mass at zero momentum or the pole mass, within the accuracy to which we 
aspire. 

Long ago, a form of the pinch technique was used to estimate the dynamical gluon mass P, 0], and interpreted 
in the factorized form outlined above. In the formulas of Q we ignore the running of the mass, and the non-linear 
integral equation then becomes: 

[g'd(q)]-^ = q'bZ - ^ I d^kH{k)H{k + q){q^ + ^} + C + . . . (12) 

where the constant C summarizes the seagull graph and other momentum-independent terms, and we also omit two- 
loop contributions. Rather than trying to solve this equation we simply replace H hy a free massive propagator, 
yielding: 

2 

[g^diq^)]-' = q'bZ + bJ{q'; Aly){q^ + + C + . . . . (13) 
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where h-uv is an ultraviolet cutoff. Omitted terms are of higher order in a dressed-loop expansion. The integral 
J{q^-Kly) is: 

J{q^;Kly) = ^ I d^k _ ^2 _^ ^^^^^^ _ ; (14) 

and the UV cutoff is defined through the Feynman-parameter representation: 

J{q^;kly)^ da\n{ ' }. (15) 

By appropriate choice of Z we make the combination Z + J finite, and define a renormalized integral J{q^\£_). Aside 
from its Feynman-parameter form J has a dispersive representation: 



The PT inverse propagator in finite terms is: 

~ 771 

[52d(g2)]-l = ^25^ + 5(^2 ^ _) J(^2. ^) + C + . . . (17) 

Note that this is consistent with the renormalization-invariance of Eq. ([T]) . We make a choice of ^ (or Z) that defines 
what we mean by A, the finite QCD scale, by requiring that 

d-'(9')- 6.9V ln(^)[l + 0(1)] (18) 

as q^ approaches infinity in any direction. (The non-leading terms are C'(lnlng^) and can affect the definition of 
A at any particular momentum, but since we deal here only with one-dressed-loop quantities we cannot use such 
higher-order terms in the analysis; one should think of A as applying to a specific range of large but finite momenta 
and that effectively incorporates terms not vanishing at infinity.) We choose: 

e - e-2A2 (19) 

and then we can set Z ~ 0. 

Since the inverse propagator is to vanish at q^ = rn^ we can eliminate C by writing: 



[q'diq')]-' = 6{J(9^;0(9^ + ^) - J(m^;e)^}. (20) 



This can be written in dispersive form by using Eq. (|16p : 



[g^d{q^)]-^^h{q^-m'){2 + \u{ — )^ daJ ^[^L + ^ ^]}. (21) 



Jim^ » a — q^ — ie a — m?) 

We now assume, as discussed in connection with Eq. (|3|), that 

g^dV) = , (22) 

which yields: 



[6.g2(q2)]-i = 2 + ln(— )- / daJ [!L + ], (23) 



A^ V fJ a — q^ — ie a \1{(7 — mn?) 

Provided that m/A is sufficiently large, a condition that we will investigate below and assume for now, the dispersion 
relation for the propagator has the properties discussed in Sec. Ill Bl above. We write: 

f{<f) - - r d-^^ (24) 
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with p a function easily read off from Eq. ((23)) : we need not record it explicitly. This equation shows that p{a) is 
positive, as we expect. An elementary calculation shows that the dispersion relation for —g^d{q^), as taken from 
Eq. dig), is: 

- g^d{q') = ^ . _ i / da- X^"^ 2 ■ ^ (25) 

m"^ — q'^ — i€ IT J (,cr — m^j((T — q^ — le) 

where i? is a positive residue, the on-shell value of (f': 

R^ifda^^. (26) 



Although each of the two factors g^iq^) and llirn^ — q^) have positive imaginary parts, the imaginary part of their 
product, which is: 

Im[-g''d{a)] = R5{a - m^) - (27) 

has one positive term from the pole and another term from the cut that is everywhere negative. This is, of course, 
required by the large-g behavior of the product which requires that the integral of the imaginary part vanish. 

One implication of the dispersive form in Eq. ([23]) is that g'^ is positive everywhere where it is real, that is, in the 
region — oo < q^ < Ann? . Another is that g^ is monotonically increasing as q^ decreases. Since we expect to be 
monotonically decreasing as the mass increases, there is a critical mass TOc such that if the physical mass exceeds rric 
there is no spurious pole, while there is such a pole if m < rric. This critical mass is determined by positivity of the 
running coupling just below the threshhold and yields mc/A — exp[\/37r/33] — 1.18. 

All the integrals in Eq. (|23p can be evaluated, giving the running charge explicitly. In the regime < < Arm? the 
explicit result is: 

[bfm-' = I, , + ^]J(?^;e) - ^J{m';m (28) 




g2 q2 2V37r 
arctan i / - — — — } 

y 4m2 - gf2 II ' 

which has an obvious analytic continuation to other regimes. 

Let us compare this result for the running charge to the old ansatz of [H, which is: 

[6g(g2)]-i=ln[(4m2-92_,e)/A2]. (29) 

It is clear that this is not accurate for the above-threshhold region q^ > Am? , because it has a pole at q^ — Am? — h? , 
which is timelike, provided that m > A/2, and lies below threshhold in the region where the running charge is real. 
However, this expression does not differ very much from the improved PT value of Eq. (j28p above in the spacelike 
regime. In Fig. [T] we compare the old expression of Eq. (|29p and the new expression in Eq. plotted vs. jm? ^ 
at a mass ratio m/A = 1. For any other value, simply add \vL(m? j to both expressions. They differ by about 15% 
or less from their average in the spacelike regime but separate increasingly for < (/^ < 477^2 ^ 1982 expression 
approaches its timelike divergence. 

For most phenomenological applications one is interested in the running charge at a small spacelike momentum 
transfer. Either the new expression from Eq. ((23|) or the old one from Eq. ([29]) shows that ^^(O) increases as m 
decreases. A single formula applies to both cases, with one parameter p whose value depends on whether the new or 
old expression is used. We find: 

Anh ln(^^) 

which is positive provided that m/A > p^^. For the new expression: 

p-,i^=exp[-l + i^(l-^)]«0.41 (31) 

and for the old fit Poj^ = 0.5 — not much different. If m/A > the squared running charge is positive for all spacelike 
(g2 < 0) momenta. 

What happens when m < m^. Is this unacceptable, or is there a fix through APT? 
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FIG. 1: Comparison of new and old expressions for [bf{q^)]-^ below threshhold. The upper curve is the new [Eq. (|23|l ] and 
the lower the old [Eq. (|29|l ] expression at m = A; the x-sais is /m? . 



III. 



APT AND OTHER MODELS 



Massless APT [T^ begins with ordinary perturbation theory for the running charge, which at one loop is: 



F{q^)^[g{q%]-^ = bH- 



-q — le 



)• 



(32) 



The tachyonic pole at q^ = —A^ is removed using a renormahzation- group- improved extension of fifty-year-old tech- 
niques that impose correct analyticity properties on certain gauge-invariant quantities such as the Adler D-function 
or the photon propagator, which amounts to postulating the dispersion relation: 



1 



I APT 



da- 



'Im F 



1 



a — q^ — le 



(33) 



As needed, the imaginary part of g\prp is positive. Since the dispersion integral can have (by fiat) no pole, the simple 
result is: 



1 



A2 



1 



47r6 4n[(-92-ie)/A2] + q- 



(34) 



with a zero-momentum value of as{0) — l/(47r&) ~ 1.4 (for three- flavor QCD). Higher-order renormalization-group 
improvement changes this value only slightly. While this value for as(0) is certainly in the right ballpark, there is an 
uncomfortable flaw in APT. It predicts that the slope das{q^) / dq^ at = is negative infinity, which is certainly 
unphysical. Consequently the predicted value in Eq. (|34p is not reliable, although that is not our main concern here. 

Obviously this massless APT treatment can be trivially extended to the old proposal of Eq. (|29p . but this just 
transfers the infinite slope to threshhold. This is a fault to be associated with the proposed running charge, which is 
simply not physical near threshhold. 

Shirkov later [ll| proposed to put in, by hand, a gluon mass in a different way that actually is close in spirit to 
the version of the PT we use here. Although Shirkov's work is unclear on some relatively insignificant details, when 
quarks are omitted it is essentially equivalent to the following expression: 



lave an unphysical pol 

with m/A > 2 this pole lies in the region A^ < q^ < Ara^. The locus of zeroes in J{q^;C) in Eq. ((1 

2{1 — 7^"^ arctan7} 



(35) 

Just as in perturbation theory, J(q^]£,) may have an unphysical pole coming from a zero of J; in the massive case 

is: 



In(^) 



(36) 
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where 



g2 



^^V4^^^^' ^^^^ 

This yields rric^ = 1. Of course, it might be possible to remove this singularity for m < rric by the same techniques 
used for massless APT, with the dispersion relation: 



hg\q') = — / da ' '"^ . (38) 

However, from the expression for J above threshhold: 

J[q^-0 = In(^) + 7"' - ^^1-' (39) 

one sees that at the critical mass both the real and the imaginary parts of J vanish at threshhold, leading to a singular 
running charge at threshhold. Presumably this is unphysical. There are no singularities for larger values of m/A, so 
it appears that for mass-improved APT there is a critical mass: mc/A = 1. 

IV. TOY VERTEX MODELS AND POSITIVITY 

It is much too difficult to consider the full Schwinger-Dyson equations even at the lowest loop level for NAGTs, so 
we construct an analytically-soluble toy model. This new model is in the spirit of the old toy model of Ref. [lH, which 
is not analytically-soluble. Both models have the same large-momentum behavior, showing asymptotic freedom and 
a beta-function with all terms negative and with factorial growth. 

Both models exploit the fact that, just as the product g^d{q^) is not only gauge-invariant but renormalization-group 
invariant, there is a similar combination for the PT vertex. Introduce the notation: 

d{q') = H{q^)Z-\q'). (40) 

With our factorization conjecture for the PT propagator this is equivalent to: 

Because g'^{q^) is positive for spacelike (negative) g^, so is Z{q^). Call the proper vertex function in the pinch technique 
r°^'^(g'i, 52, ga)- This vertex function, which like the PT propagator is gauge-invariant and process- independent, obeys 
a Ward identity of QED type, with no contributions from ghosts, schematically of the form qi - T — Z{q2) — Z{q3). 
The gauge- invariant and renormalization-group invariant we call G (irrelevant group and spin indices omitted): 

/V/ ^ 31^(91,92,93) 

G(qi,q2,q3) — — ^- (42) 

{Z{q,)Z{q,)Ziq,)y/^ 

where T is the PT proper vertex function (again, irrelevant indices omitted) and Z is a factor in the propagator, as 
given in Eqs. (|40I41[) . When all the momenta are 0{q) the Ward identity tells us that F ~ 9^g~^ at large momentum, 
and then Eq. ()42p shows that G ~ g{q), as would be expected for a gauge-invariant vertex function. 

One might think that the PT Schwinger-Dyson equation (SDE) for T explicitly involves the PT propagator d that 
has a non-positive imaginary part. Instead we remark that this equation can be rewritten in terms only of the normal 
propagators H and the special vertex G. This is important because d{q'^) itself violates the K-L representation, and 
if the skeleton graphs of the SDE were to be modeled by replacing bare propagators by d there could possibly be 
positivity problems in the SDE arising from the d terms. Schematically the one-dressed-loop SDE is: 

G = Go + J G'^H^ + ... (43) 

Here G is the Born term, behaving like (Ing^)"'^/^ when all momenta are large and 0{q). Note that this equation is 
independent of the coupling constant g, as it must be if G is renormalization- invariant; this independence holds for 
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all vertex skeleton graphs. We can now draw conclusions based on the (Euclidean) positivity of the H propagators 
without fear of difhculties arising from non-positivity of the spectral function for the propagator itself. 

So far we have not considered numerator factors. Neither toy model has them, but each roughly accounts for them 
by dropping one of the propagator factors in Eq. We know that G ^ (In 5^)"-'^/^ at large q, which dominates over 
the inhomogeneous Born term. Consequently, in the toy model this inhomogeneous term is dropped. Furthermore in 
the toy model, G depends on only one momentum, and only the one-loop skeleton graph is saved. The original toy 
model equation is then: 

^^'^^ ^2^1 (fc2 - m2 + ie)[{q + kf - m2 + ie] ^'^''^ ^^^^ 

where h is the usual (no-quark) one-loop coefficient in the beta-function. This is a universal equation for any coefficient 
h in the beta-function, as one sees by using the vertex R = b^/'^G in place of G. 

We will work in Euclidean space, defining = as the Euclidean square of the momentum, positive for spacelike 
q. Then the kernel of the non-linear integral equation (j44p is positive. If m 7^ the kernel is nowhere singular, but if 
•m — Q the kernel is singular at zero momentum. One implication is that the massless G is necessarily zero or singular 
at zero momentum, already suggesting the necessity of a mass. We briefly review these facts for the massless model. 



A. Massless toy model 

The massless toy model can be converted [l^ to a differential equation: 

Gu + Gt = -lG'' (45) 

where the subscripts indicate derivatives with respect to the variable t = \n{Q^ /A^) and A is the usual QCD mass 
scale. We would like to impose physically-sensible boundary conditions at — 0, or t — —00, but this is impossible: 
The massless version of Eq. (H^ is singular at zero momentum unless G(0) = 0, which we forbid. Nevertheless the 
differential equation can be solved, showing features expected from perturbation theory, and we can impose boundary 
conditions aX t > and study the ultraviolet behavior. 

In the ultraviolet regimes one finds results familiar from the RG: The functional form of the asymptotic vertex is 
precisely that of the full NAGT except for the value of some numerical coefficients, and all signs agree with what is 
needed for asymptotic freedom. For large t the second derivative term is non-leading; if dropped, the general solution 
to the first-order differential equation is: 

bt = const. (46) 



This coupling is singular, as massless perturbation theory must be. When the second derivative term is kept, a 
solution is generated which has all the same terms as the all-order perturbative running charge in true QCD, but with 
somewhat different coefficients. All these coefficients have the correct sign for an asymptotically-free theory. 

There is also an interesting beta-function, governed by its own differential equation. This comes from the relation: 



plus the equation (|^5)) for the vertex, and is: 



Gt = ^/3(G) (47) 



+ = (48) 



It was shown [1^ that this beta-function behaves like —g^N\{bcg^)^ for some positive constant c, qualitatively just 
the same as in any asymptotically- free NAGT, and that the beta- function solving Eq. (|48)) is singular at a finite upper 
critical coupling Qc- 



B. Massive toy model 



We give here a new toy model inspired by but differing slightly from the original massive toy model of Eq. (j44p . 
The new model is probably about as accurate as the original in modeling the true SDE, but it can be analytically 
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solved and shows a critical mass value. It has the same ultraviolet (massless) properties as the old model and as QCD 
does itself. As before we take as a free massive propagator. There is no longer a simple differential form of the 
original toy model equation Eq. (|44p when masses are included, but the following approximation to the Euclidean 
angular integration does give an ordinary differential equation: 



27r2 [(Q + X)2+to2] Q2 + m^ 

The approximation is exact for large Q,K as well as when Q > 0, K — (or K > 0,Q — 0) and is otherwise too 
large by a factor which is at most about 1.3 times the true angular integral a.t Q — K — m. From this follows the 
one-dimensional integral equation: 

G(q2) ^ Ji(g2) + J2(g2) (50) 



h r 



2{Q^ + m^)Jo if2 



(51) 



Every solution of this integral equation G satisfies the differential equation: 

[(Q2 + m')G{Q')r' ^ - ^SfKi (52) 



where a prime indicates differentiation with respect to . If we now define t as: 



, , , Q' + m2 

t = H — ), (53) 

the differential equation can be written as an extension of the massless equation: 

Gu + G, = -\gHI-'^). (54) 

The difference from the massless equation is that at Q = the variable t is finite, not — oo. 

It is not difficult to check from the integral equation (|50p that at large momenta only the J2 term is leading, 
yielding the exact leading behavior G(g2) ^ l/(61n02)i/2 at large Q^. The Ji term in Eq. (HDD is 0[(lnQ2)-3/2)] 
and non-leading. In neither the original toy model nor at present are we interested in such non- leading terms, so we 
define our new toy model by dropping the Ji term. Note also that Ji = 0{Q'^) at small , while J2 = 0(1), since 
we require that (5(0) ^ 0; the missing Ji term is non-leading in the infrared as well. [One can show that the 0{Q'^) 
terms in Ji and J2 cancel in the full equation, leaving corrections to (5(0) of 0{Q^). In fact, the exact solution G of 
Eq. (j50[) with both terms show the self-consistent behavior leading correction at small momentum: 

(?(Q2)^(?(0)-^(^)3 + ...] (55) 
12 

Evidently the resulting integral equation with no Ji term satisfies a first-order differential equation, which is just 
Eq. (|5i|) without the Gtt term. This equation has the exact solution Gi. 

= 72 ) + , ^2 ]- (56) 



(72(Q2) ^ ^ A2 ^ Q2+m2J- 

One can check that the beta-function coming from Gi has not only the usual —hG\ term but also terms involving 
non-perturbative quantities such as exp[— 1/(6G2)] This approximation already shows a mass mo at which as(0) is 
singular: toq/A = e~^/2 ~ 0.61; smaller values lead to a pole in Gi. (We have done a quick numerical study of the 
differential equation Eq. ([51]) of the modified model, and found too /A « 0.66.) The actual critical mass nic might be 
about twice as large as mo, based on our experience with the propagator model. 
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-0.5- 



-1.5- 
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FIG. 2: Solutions to the differential equation (|57|l for K—1 (right-most curve) to K=5 (left-most curve). Here x — Vbg and 



TABLE L Upper-bound couplings for various K 



K = 


1 


2 


3 


4 


5 




3.1 


1.3 


0.90 


0.70 


0.58 



C. Is a gluon mass really necessary? 



One may ask whether higher-order effects in a massless theory can somehow do away with the need for a gluon mass. 
To study this possibility, let us now extend the basic one-loop beta-function equation (|48p to mimic multi-loop vertex 
contributions. Of course, we can only aspire to qualitative accuracy, looking for mechanisms of limits on rather 
than for accurate values of these limits. These models are essentially those of [40] devised for (f>g, with coupling G. In 
this theory it has been shown that the sum of all iiT-loop vertex graphs is always positive and grows like Kl{cG^)^ , 
where c is a positive number. We conjecture that the same holds for d — 4 NAGTs, and that an appropriate extension 
of the beta-function differential equation to higher loops is: 



1 



K 

^ j=i 



where 



= y(x) = \/6/3(g) 



(57) 



(58) 



and the prime denotes an x derivative. One can verify that the beta-function solving this equation grows factorially 
like Kl{cg^)^ for some c, whose value is irrelevant to our investigation. Our only interest is in seeing what happens 
when the right-hand side of the vertex Schwinger-Dyson equation has terms that grow factorially with loop number 
K. In fact, we have tested the sensitivity of Eq. ((57)) to fairly major changes in the coefficients on the right-hand side, 
and find little change from the solutions to Eq. ([57]) as shown in Fig. [21 

Table 1 show the upper bounds ac for quarkless QCD, defined as the value of gl/{4:Tr) at which the beta-function 
crosses the real axis, for various values of K. One should expect the ac values to decrease as K increases, since 
otherwise the right-hand side of the differential equation will grow too large to balance the left-hand side. For all 
values of K there is a singularity of the form /3 ^ {g^ — g)^^"^, where f3 has infinite slope at g = gc and then turns 
imaginary for g > gc- Note that there seems to be numerical convergence toward a value near 0.5, not far from "best" 
estimates based on the one-dressed-loop pinch technique with a mass. But there is always a singularity. 
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Perhaps some way of summing the non-Borel-summable series of Eq. (|57p would remove this singularity. We have 
tried, again in the spirit of "regulating" the all-orders behavior with a principal part integral form of the Borcl 
integral corresponding to the sum in Eq. ([ST]) , using: 

y{l + ^y')=L{x). (59) 

with: 

^(^) = -472 / — — }• (60) 

' Jo a ' I — a 

The equation for the vertex itself, analogous to Eq. (I45p . is: 

Gu + Gt^b-'/^Lib'/^G). (61) 

Actually, some form of principal-part regulation is demanded by dynamical boson mass generation [40| . but we need 
not inquire further into that here. We can understand the basic behavior of the equation by looking at the degree to 
which factorial growth at large coupling is tamed by the principal-part prescription. The power-series expansion of 
L{x) is 

L{x) = - ^ E ^'""^'m -h+ 0(e-''^') (62) 

^''^ K=l 

and its asymptotic behavior at infinity is —2x. This is very different from the finite- -fC models of Eq. (|57p . In these, the 
increasingly-strong growth with x of the right-hand side as K gets larger means that any singularity occur at smaller 
values of x. But for equations with the L{x) source there is no such movement toward smaller couplings because 
L{x) is not growing rapidly at large a;. In the beta-function equation (j59p one easily finds the large- a; asymptotic 
behavior y{x) — > — ca; with c = — 1 + \/5- So instead of generating a finite-g singularity, the beta- function turns 
from —bg^ behavior near the origin to linear at large enough g; numerical simulations confirm this. For the massless 
vertex equation (|6ip there is still a singularity in the infrared, so that G gets large and the large x behavior of L{x) 
matters. In this regime the massless equation (I6ip becomes linear, and one finds unphysical behavior of the type 
(ln(5^)~^/^cos(\/31n(5^). So massless vertex dynamics is not regularized by the specific behavior of L{x). Mass is 
important not only for the right-hand sides of the vertex and beta-function equations, as summarized by the function 
L, but it is important in the left-hand side of such equations, as we show below. In fact, the mass damping is so 
strong that it is probably unnecessary to worry about vertex graphs of very high order; G does not get large enough 
to probe the asymptotic Hmit of L{b&), as it does for the massless case. 



V. SUMMARY AND CONCLUSIONS 



We argue that although the gauge-invariant PT propagator does not obey the K-L representation, it is the product 
of two factors that do have the required positive spectral functions. We show that this holds true for an analytic 
approximation to the one-dressed-loop PT propagator equation, and that this approximation implies a critical mass 
such that the true dynamical mass m must exceed nic or spurious singularities arise. We construct an analytically- 
soluble toy model of one-dressed-loop PT three-gluon vertex model with a cubic non-linearity and show how the 
positivity of the factor H in the PT propagator plays an essential role in the vertex dynamics, in particular the 
occurrence of a mass value which must be exceeded by m to avoid unwanted singularities. We argue that higher- 
order, even regulated all-order, extensions of the massless toy model equations do not remove these singularities. 
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